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An  attractive  option  for  constructing  thermoelectric  generators  (TEGs)  is  to  incorporate  a  water-fed  heat 
exchanger  with  commercially  available  thermoelectric  modules.  In  this  paper,  two  different  thermoelec¬ 
tric  models  are  applied  to  predict  the  energy  conversion  performance  of  the  TEGs.  The  first  model 
employs  a  derivation  of  the  Carnot  efficiency.  The  second  model  presents  a  rigorous  interfacial  energy 
balance  by  capturing  Joule  heating,  Seebeck,  Peltier  and  Thomson  effects,  yielding  better  predictions  of 
the  conversion  capability.  This  model  is  then  used  to  perform  a  computational  examination  of  the  TEGs 
embedded  in  30  different  configurations,  which  allows  the  identification  and  quantification  of  key  design 
parameters  including  flow  types,  hot  stream  inlet  temperatures  ( TinJl0t ),  pressure  drops  (AP),  cross- 
sectional  area  (Ac),  channel  length  ( Lch )  and  number  of  channels.  The  positive  effects  of  TinJlot  and  AP 
can  be  easily  captured  and  parallel  flows  of  the  hot  and  cold  streams  are  found  to  provide  greater  overall 
TEG  efficiency  as  compared  to  counter  flows.  In  general,  micro-sized  Ac  reduces  temperature  gradients 
across  the  channels,  providing  a  greater  AT  across  the  thermoelectric  material.  However,  enhancements 
of  the  conversion  capability  are  eventually  limited  by  the  reduced  convective  heat  transport  due  to 
increased  flow  resistance.  Finally,  further  improvements  in  the  power  generation  are  achieved  by  reduc¬ 
ing  Lch  while  increasing  the  number  of  channels.  The  resulting  reduction  in  flow  resistance  is  found  to 
facilitate  increases  in  convective  heat  transfer,  as  well  as  in  AT,  and  thus  a  great  increase  in  conversion 
efficiency  ( rj ). 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectrics  are  in  principle,  attractive  for  the  conversion  of 
waste  heat  into  electric  power.  The  underlying  mechanisms  dictat¬ 
ing  thermoelectric  effects  were  not  well  understood  until  the  dis¬ 
covery  of  the  electron  at  the  end  of  the  19th  century.  It  is  now 
known  that  electrons,  or  holes,  in  a  solid  carry  both  charge  and 
heat.  When  a  temperature  drop  is  applied  across  a  thermoelectric 
material,  as  shown  in  Fig.  1,  this  will  cause  the  net  diffusion  of  elec¬ 
trons  (holes)  from  the  hot  side  to  the  cold  side.  The  electrical 
potential  produced  by  the  temperature  difference  is  known  as 
the  Seebeck  effect,  and  devices  that  generate  electrical  power  by 
virtue  of  the  Seebeck  effect  are  called  thermoelectric  generators 
(TEGs).  Alternate  electrical  connections  of  the  hot  sides  and  cold 
sides  of  a  series  of  n-  and  p-type  materials  allow  the  build-up  of 
any  desired  voltage. 
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TEGs  have  great  potential  in  the  direct  conversion  of  waste-heat 
energy  from  power  plants  and  automotive  vehicles  into  electric 
power  where  it  is  unnecessary  to  consider  the  cost  of  the  thermal 
energy  input,  and  offer  reliable  power  in  remote  areas  such  as  in 
space  and  at  mountain  top  telecommunication  sites  [1  ].  Unlike  tra¬ 
ditional  dynamic  heat  engines,  TEGs  present  several  distinct 
advantages  over  other  technologies,  such  as  being  simple,  compact 
and  highly  reliable,  being  capable  of  operating  at  elevated  temper¬ 
atures,  being  suited  for  small-scale  and  remote  applications,  and 
being  environmentally  friendly  [2,3].  A  major  limitation  is,  how¬ 
ever,  the  low  conversion  efficiency  of  current  thermoelectric 
devices  (typically  around  5-10%  [4]),  which  is  mainly  caused  by 
the  available  thermoelectric  materials  in  the  direction  of  figure  of 
merit  and  has  restricted  their  use  in  specialized  fields  with  exten¬ 
sive  applications.  Since  the  promotion  of  the  intrinsic  efficiencies 
of  the  TE  materials  at  device  level  has  yet  to  be  proved,  there  are 
ongoing  attempts  to  increase  the  competitiveness  of  TEGs  by 
improving  the  way  in  which  they  are  currently  used,  offering  some 
significant  guidelines  for  design  optimization  [5-7].  The  maximum 
efficiency  of  TEGs  is  subject  to  the  constraints  of  the  second  law  of 
thermodynamics,  with  the  limit  being  expressed  as  r\max  =  A T/T^, 
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Nomenclature 

X 

axial  distance  along  fluid  flow  axis,  m 

G 

electric  conductivity  of  TE  material,  S/m 

L 

TEG  length,  m 

5 

Seebeck  coefficient  of  TE  material,  V/I< 

W 

TEG  width,  m 

P 

density,  kg/m3 

H 

TEG  height,  m 

P 

pressure  in  flows,  Pa 

Subscripts 

T 

temperatures,  K 

h,  c 

hot  side  and  cold  side  of  thermoelectric  layer 

cP 

specific  heat  capacity  at  constant  pressure,  J/kg  K 

TE 

thermoelectric 

Nu 

Nusselt  number 

ch 

channel 

A 

area,  m2 

avg 

average 

G 

volume  flow  rate,  m3/s 

in_hot, 

in_cold  inlet  of  hot  and  cold  stream 

1-5 

length  number 

Greek  symbols 

conv 

convection 

K 

thermal  conductivity,  W/K  m 

Heat  absorbed 


Fig.  1.  Schematic  illustration  of  a  typical  TEG.  The  construction  consists  of  a  pair  of 
p-type  and  n-type  semiconductor  materials  forming  a  thermocouple.  Electron 
(Hole)  diffusion  is  shown  in  this  figure  driven  by  a  temperature  difference. 

where  Th  is  the  temperature  of  the  thermoelectric  module  hot  side, 
and  a  principal  focus  has  been  on  improving  generator  efficiency 
by  maximizing  the  temperature  drop  (AT)  across  the  thermoelec¬ 
tric  modules.  Traditional  efforts  have  aimed  to  achieve  this 
through  the  design  of  heat  exchangers.  Esarte  et  al.  [8]  theoreti¬ 
cally  addressed  the  optimization  of  heat  exchanger  geometries 
and  various  operating  conditions  including  fluid  flow  rates,  fluid 
properties  and  inlet  temperatures.  Crane  et  al.  [9]  investigated  a 
computational  model  for  the  integration  of  a  thermoelectric  heat 
exchanger  in  a  cross  flow  configuration.  In  the  case  of  TEGs  with 
large-scale  multi-panels  [10]  and  cylindrical  multi-tubes  [11],  a 
mathematical  expression  for  power  output  was  deduced,  consider¬ 
ing  different  geometrical  systems.  Experiments  have  also  been 
conducted  on  TEGs  intended  for  low-temperature  waste  heat 
power  recovery  to  assess  effects  of  inlet  temperature  and  flow  rate 
on  power  output  and  conversion  efficiency  [12]. 

All  of  these  research  efforts  combined  thermoelectric  modules 
with  macro-scale  heat  exchangers,  but  growing  applications  such 
as  autonomous  micro-systems  or  wearable  electronics  look  for 
micro-scale  power  generators.  A  micro  TEG  was  proposed  to  pro¬ 
vide  electric  power  for  an  electronic  chip  in  a  domestic  gas-moni¬ 
toring  system.  In  excess  of  1.5  V  could  be  produced  when  AT  of  a 
few  tens  of  degrees  was  established  [13].  Moreover,  the  utilization 
of  TEGs  in  reciprocating  internal  combustion  engines  is  another 


novel  application.  A  comprehensive  theoretical  study  of  applying 
a  TEG  to  automotive  engine  waste  heat  recovery  is  conducted 
[14].  By  utilizing  AT  of  563  K  between  a  micro-channel  heat  sink 
and  exhaust  gas,  the  maximum  power  of  51.13  mW/cm2  can  be 
reached. 

Enormous  amounts  of  research  have  been  conducted  on  appli¬ 
cations  of  TEGs.  However,  these  thermoelectric  systems  are  con¬ 
strained  by  their  low  thermodynamic  efficiencies,  which  means 
that  a  comparatively  large  amount  of  heat  is  required  to  produce 
a  given  quantity  of  electricity.  Thus,  there  is  a  clear  need  for  further 
improvements  in  the  design  of  the  heat  exchangers  as  well  as 
matching  of  the  heat  exchanger  and  thermoelectric  designs.  In  this 
paper,  the  device  performance  of  a  water-fed  heat  exchanger  based 
TEG  in  a  range  of  scales  is  estimated  by  a  simulation  approach, 
which  offers  more  detailed  predictions  than  analytical  methods, 
and  can  be  more  cost-effective  than,  or  eventually  help  guide, 
experiments.  The  simulations  are  performed  using  the  commercial 
FEA  package  COMSOL™  [15].  Two  different  thermoelectric  models 
are  applied  in  this  paper.  The  first  is  the  Simple  model  derived  from 
the  Carnot  efficiency,  and  the  second  is  the  Coupled-field  model 
presenting  a  rigorous  interfacial  energy  balance  by  capturing  Joule 
heating,  Seebeck,  Peltier,  and  Thomson  effects. 

Experiments  conducted  by  Niu  et  al.  [16]  are  employed  here  to 
demonstrate  the  applicability  of  the  two  thermoelectric  models  for 
evaluating  TEGs  device  performance.  The  model  showing  the  best 
agreement,  by  comparison  with  the  experimental  measurements, 
is  used  to  perform  a  computational  examination  of  TEGs  embedded 
in  30  different  configurations,  which  allows  the  identification  and 
quantification  of  systematic  effects  of  key  design  parameters 
including  flow  types  (parallel  flows  &  counter  flows),  hot  stream 
inlet  temperatures  ( Tin  hot ),  pressure  drops  (AP),  cross-sectional 
area  (Ac),  channel  length  ( Lch )  and  number  of  channels.  The  positive 
effects  of  Tm  hot  and  AP  on  energy  extraction  are  captured.  Micro¬ 
scale  Ac  values  provide  smaller  temperature  gradients  across  the 
channels  and  thus  a  greater  AT  across  the  thermoelectric  material. 
However,  enhancements  of  the  conversion  capability  are  eventu¬ 
ally  limited  by  reduced  convective  heat  transport  due  to  increased 
flow  resistance.  Finally,  further  improvements  in  the  power  gener¬ 
ation  are  achieved  by  reducing  Lch  and  increasing  the  number  of 
channels. 

2.  Thermoelectric  models 

2.1.  Simple  model 

In  the  Simple  model  established  previously  by  the  authors  [17], 
the  three-dimensional  Navier-Stokes  and  energy  equations  com¬ 
bined  with  the  continuity  equation  are  solved  simultaneously  to 
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determine  the  temperature  fields  across  the  thermoelectric 
module. 

The  governing  equations  read  as  follows, 

Continuity  equation: 


Momentum  conservation  equation  reads, 
p(u  •  Vu)  =  -VP  +  V  •  cp  +  F  (9) 

Energy  conservation  for  the  liquid  domain  shows, 


V  •  (pu)  =  0 

Momentum  conservation  equation: 
p(u  •  Vu)  =  -  VP  +  V  •  cp  +  F 
Energy  equation: 


u  VT  = 


pCp 


v2t 


_Q_. 

pCp 


s 

pCp 


(1) 

(2) 

(3) 


where  F  represents  body  forces  acting  on  the  fluid,  cp  is  the  stress 
tensor,  5  is  the  dissipation  function  due  to  the  viscous  force,  and 
Q  is  the  rate  of  internal  heat  generation  within  the  solid  domain. 
The  energy  equation  is  solved  in  the  fluid  and  solid  domains  where 
the  heat  transfer  is  strictly  by  convection  and  conduction  respec¬ 
tively.  For  the  fluid  domain,  Q  is  zero;  for  the  solid  domain,  the  con¬ 
vective  term  and  5  are  zero,  while  Q  =  -  W' .  W'  is  the  electric  power 
extracted  by  the  TEGs  and  can  be  expressed  as, 


dW'  =  qm3Xdq 


(4) 


q  is  the  heat  flow  rate  by  conduction  across  the  thermoelectric 
layer,  which  may  be  evaluated  from  Fourier’s  law, 


dq  =  - KTEdydz ^  -  KTEdxdy ||  -  KTEdxdz || 


(5) 


Limited  by  the  second-law  of  thermodynamics,  the  ideal  effi¬ 
ciency  of  a  TEG  operating  as  a  reversible  heat  engine  is  Carnot  effi¬ 
ciency,  given  by 


^Carnot  ^  y 

1  h 


(6) 


The  maximum  conversion  efficiency  of  an  irreversible  TEG  can 
be  estimated  by  [18] 


max  ^7  Carnot 


1+Z^-l 


yi+zi^+t 


(7) 


u-VT  =  Tv2T  +  T  (10) 

pCp  pCp 

Energy  conservation  for  the  solid  domain  involves  Joule  heating, 

v(kte-vt)  — r-j  -(J0-r  +  /oTE-j  =  o  (ii) 

In  this  model,  the  equations  of  heat  flow  Eq.  (11)  and  of  conti¬ 
nuity  of  electric  charge  [19]  Eq.  (12)  are  coupled  by  the  set  of  ther¬ 
moelectric  constitutive  equations,  Eqs.  (13)  and  (14). 


VJ  =  0  (12) 

q  =  [n]  -  J-  [k]  •  vr  (13) 

J  =  [°1 '  (E  -  [<5]  •  VT)  (14) 


where,  q  is  the  heat  flux  vector,  J  is  the  electric  current  density  vec¬ 
tor,  E  is  the  electric  field  intensity,  [k]  is  the  thermal  conductivity 
matrix,  [a]  is  the  electrical  conductivity  matrix,  [S]  is  the  Seebeck 
coefficient  matrix,  [IT]  =  T[S]  is  the  Peltier  coefficient  matrix.  In 
the  absence  of  time-varying  magnetic  fields,  E  is  irrotational,  and 
can  be  derived  from  an  electric  scalar  potential  cp , 

E  =  -\7cp  (15) 

The  flow  field  needs  to  be  incorporated  into  the  thermoelectric 
schemes,  which  requires  the  addition  of  Navier- Stokes  model.  The 
equations  governing  weakly  compressible  Navier-Stokes  Flow  can 
be  found  in  Eqs.  (l)-(3). 

Therefore,  the  electric  power  W'  can  be  calculated  as 


W  = 


vl 


2(PTe  +  Ric 


(16) 


where,  V0c  is  the  open  circuit  voltage,  PTe  is  the  internal  thermo¬ 
electric  resistance,  and  Rioad  is  the  electric  resistance  of  the  external 
load.  For  the  maximum  electrical  power,  R[oad  =  RTE. 


2.3.  Conversion  efficiency 


The  figure  of  merit  Z  for  the  thermoelectric  devices  can  be 
defined  as, 


In  general,  a  TEG  exhibits  low  efficiency  due  to  the  relatively 
small,  dimensionless  figure-of-merit  (ZT,  a  metric  of  measuring 
currently  available  thermoelectric  materials’  performance  at  a  par¬ 
ticular  temperature).  By  solving  Eqs.  (l)-(8)  simultaneously,  the 
temperature  profiles  of  thermoelectric  layers  are  first  obtained 
and  substituted  into  Eq.  (4)  to  get  the  value  of  W'.  This  value  is  then 
coupled  back  to  Eqs.  (l)-(8)  to  get  the  updated  temperature  field, 
and  the  power  W'. 

2.2.  Coupled-field  model 

While  the  Simple  model  has  being  frequently  used  to  estimate 
the  potential  performance  of  thermoelectric  devices,  the  diversity 
and  complexity  of  thermoelectric  applications  necessitate  a  fully 
Coupled-field  model,  which,  in  addition  to  Joule  heating,  accounts 
for  Seebeck,  Peltier  and  Thomson  effects  as  coupling  mechanisms 
between  thermal  and  electric  fields.  The  equations  governing  the 
multidimensional  temperature  and  electrical  potential  distribu¬ 
tions  in  the  TEGs  under  steady-state  conditions  and  in  the  absence 
of  an  applied  magnetic  field  are  shown  as, 


A  term  r\  is  referred  to  as  the  overall  TEG  efficiency  and  can  be 
defined  as, 


W  _  W 


(17) 


3.  Solution  procedure 

3.1.  Model  evaluation 

While  the  simulation  is  a  common  engineering  practice,  the 
validity  of  the  proposed  methodology  becomes  very  important. 
The  experiments  conducted  by  Niu  et  al.  [16]  are  used  to  demon¬ 
strate  the  applicability  of  each  thermoelectric  model  for  conversion 
capability  prediction.  The  variations  of  power  generation  with 
Tin_hot  from  both  simulations  and  experiments  are  shown  in 
Fig.  2.  The  maximum  power  output  obtained  using  the  Simple 
model  is  approximately  30%  higher  than  that  of  the  experiment, 
which  may  due  to  the  electric  field  is  not  considered  in  this  model. 
The  comparison  between  model  predictions  and  experiments 
clearly  indicates  that  the  Coupled-field  model  shows  better  perfor¬ 
mance  as  compared  to  the  Simple  model,  which  is  then  used  in  the 
succeeding  simulations  to  examine  the  key  design  parameters  for 
TEGs. 
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Fig.  2.  Model  validations  with  maximum  power  output  at  the  reference  conditions  (rin_co/d  =  293  K,  Gin_hot  =  0.4  m3/h,  Gin_co/d  =  0-3  m3/h). 


L,  30mm 


Pressure  Inlet,  Hot  stream 


Pressure  Inlet,  Cold  stream 


Fig.  3.  Schematic  layout  of  the  proposed  TEG  considered 


in 


this 


study 


H.  14mm 


with  boundary  specifications. 


Outlet 


3.2.  Simulation  specification 

In  this  study,  considering  a  basic  TEG,  an  entire  thermoelectric 
system  of  width,  height  and  length  12  mm  by  H  14  mm  by  L 
30  mm  is  modeled  consisting  of  a  water-fed  heat  exchanger  and 
commercially  available  thermoelectric  modules.  The  numerical 
heat  exchanger  model  is  integrated  with  models  for  thermoelectric 
modules,  which  are  positioned  between  hot  and  cold  flow  streams 
as  shown  in  the  schematic  of  a  TEG  section  in  Fig.  3.  The  heat 
applied  from  a  hot  stream  is  injected  at  the  top  of  a  thermoelectric 
layer  and  extracted  by  a  cold  stream  at  the  bottom.  The  Seebeck 
effect  can  therefore  generate  electric  power  from  AT  across  the 
thermoelectric  layer.  All  subsequent  thermoelectric  schemes  here 
are  analyzed  based  on  the  Coupled-field  model  at  a  steady  state 
operation,  utilizing  the  commercial  FEA  package  COMSOL™ 
through  the  Physics  Interface  Builder.  The  fluid  flow  in  this  study 
is  along  the  x-axis  shown  in  Fig.  3,  and  is  assumed  to  be  incom¬ 
pressible.  The  continuum  assumption  is  accepted  for  the  fluid  stud¬ 
ied  here,  as  it  is  liquid  water.  The  no-slip  and  no  flow-through 
boundary  conditions  on  all  fluid/solid  interfaces  are  invoked.  Hot 
water  enters  inlet  ports  at  373  K,  while  cold  water  enters  at 
293  K.  Wall  thicknesses  are  neglected.  There  is  natural  air  convec¬ 
tion  from  all  external  surfaces,  and  the  rate  of  heat  transferred  to 
the  surrounding  by  convection  is  calculated  by, 

Q  —  ^Xconv  '  AConv  '  C^s  —  ^oo)  ("18) 


where  hconv  is  the  convective  heat  transfer  coefficient,  Aconv  is  the 
area  where  convective  heat  transfer  occurs,  Ts  is  the  surface  temper¬ 
ature,  is  the  ambient  temperature. 

Different  numbers  of  grid  cells  are  used  to  perform  a  grid  sensi¬ 
tivity  study  at  the  reference  conditions  (Tin  cold  =  293  K,  TinJl0t  = 
373  K,  and  AP  =  2  bar).  As  shown  in  Fig.  4,  a  hexahedral  mesh  with 
at  least  a  Degree  of  Freedom  (DOF)  of  700,000  is  necessary  to  ob¬ 
tain  a  grid-independent  result. 

3.3.  Thermoelectric  material  characterization 

The  thermoelectric  material  is  a  fundamental  component  of  a 
TEG.  Recent  advances  in  materials  and  materials  processing  as 
summarized  in  recent  reviews  [20,21]  have  led  to  higher  ZT  values, 
and  thus,  higher  theoretical  rjmax .  For  AT  >  600  K,  new  thermoelec¬ 
tric  materials  i.e.  Zn4Sb3-based  compounds  have  exhibited  rjmax 
near  15%  or  more  for  AT  >  800  I<  [22].  However,  Bi2Te3-based  com¬ 
pounds  offer  the  best  performance  for  Th  in  a  range  of  350-400  K 
[23]  where  inexpensive  fluids  such  as  water  can  be  used.  These 
materials  are  taken  as  the  focus  of  the  current  study,  but  our  con¬ 
clusions  are  more  generally  applicable.  The  thermoelectric  leg  in 
the  simulations  shown  in  Fig.  3  is  1  mm  by  1  mm  by  1.4  mm, 
capped  by  thin  copper  electrodes  with  the  height  of  0.3  mm.  The 
material  properties  are  tabulated  in  Table  1.  Usually  these  are  tem¬ 
perature-dependent  and  may  be  anisotropic,  but  we  take  them  to 
be  isotropic  and  constant. 
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Degree  of  Freedom  (DOF) 


Fig.  4.  Power  generation  for  the  case  depicted  in  Fig.  3  vs.  grid  size  using  the 
Coupled-field  model  at  the  reference  conditions  ( Tincoid  =  293  K,  Tin  hot  =  373  K,  and 
AP  =  2  bar). 


Table  1 

Material  properties  for  TE  modules. 


Symbol 

BhTe3 

Electrode  (Copper) 

Seebeck  coefficient 

V/K 

p:  2.2  x  1(T4 
n:  -2.2  x  10“4 

6.5  x  10“6 

Thermal  conductivity 

k,  W/(K  m) 

1.5 

350 

Electric  conductivity 

g,  S/m 

1.0  x  105 

5.9  x  108 

Density 

p,  kg/m3 

7740 

8920 

Heat  capacity 

Cp,  J/(kg  I<) 

154.4 

385 

4.  Results  and  discussion 

For  the  first  calculation,  the  temperature  profile  of  the  baseline 
TEG  depicted  in  Fig.  3  is  shown  in  Fig.  5(a),  while  the  voltage  profile 
generated  based  on  the  Coupled-field  model  is  indicated  in 
Fig.  5(b).  An  electrical  power  of  0.65  W  is  obtained.  As  might  be  ex¬ 
pected,  an  optimal  heat  exchanger  for  power  generation  would 
maximize  AT  across  the  thermoelectric  layer  by  having  minimal 
heat  transfer  resistance  and,  thus,  minimum  temperature  drops 
at  both  the  hot  and  cold  side  solid-liquid  interfaces,  as  well  as  min¬ 
imum  temperature  gradients  across  the  height  of  the  channels.  In 


Fig.  6.  Variation  of  the  power  generation  and  conversion  efficiency  with  TinJl0t 
( Tin_Coid  =  293  K  and  AP=2  bar). 


the  following  sections,  the  key  design  criteria  are  carefully 
examined. 

4.1.  Effect  of  hot  stream  inlet  temperature  (Tinhot) 

Fig.  6  shows  the  increase  in  power  generation  and  conversion 
efficiency  with  increasing  the  inlet  temperature  of  the  hot  liquid, 
Tin  hot ,  assuming  everything  else  (notably  the  cold  inlet  tempera¬ 
ture)  is  kept  constant.  An  increase  in  Tin _hot  has  two  counteracting 
effects  as  far  as  r\  is  concerned  (according  to  Eq.  (17)),  but  the 
increase  in  AT  dominates  that  and  the  overall  power  generation 
increases  faster  than  linearly.  Not  surprisingly,  it  is  thus  desirable 
to  have  the  hot  water  delivered  at  as  high  a  temperature  as  possi¬ 
ble.  Lowering  the  temperature  of  the  cold  liquid  would  have  an 
even  stronger  effect  on  r\ ,  but  that  is  not  as  often  practical. 

4.2.  Effect  of  flow  type  and  pressure  drop  (AP) 

Simulations  are  carried  out  for  both  parallel  flow  (hot  and  cold 
flows  in  the  same  direction)  and  counter  flow  (hot  and  cold  flows 
in  the  opposite  direction)  within  the  heat  exchanger.  Fig.  7(a) 
compares  the  variation  of  AT  with  flow  distance  from  the  inlet 
plane  for  parallel  flow  and  counter  flow,  with  a  channel  length  of 
2  m.  For  given  inlet  temperatures  of  each  fluid,  there  is  a  significant 


T,  K  V,  V 


(a)  (b) 


Fig.  5.  (a)  Temperature  profile  of  the  water-fed  TEG  at  the  reference  conditions  ( Tin_co/d  =  293  K,  Tinhot  =  373  K,  and  AP  =  2  bar);  (b)  Voltage  profile  generated  from  the  water- 
fed  TEG  based  on  the  Coupled-field  model. 
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Fig.  7.  Variation  of  (a)  the  temperature  drop  with  flow  distance  from  the  inlet  plane  ( Lch  =  2  m);  (b)  the  average  temperature  drop  between  two  fluids  at  the  channel  ends 
with  Lch;  (c)  power  generation  with  Lch  for  parallel  flow  vs.  counter  flow  at  the  reference  conditions  ( Tin  coid  =  293  K,  Tin  hot  =  373  K,  and  AP  =  2  bar). 


variation  in  AT  for  parallel  flow  but  a  relatively  smaller  variation 
for  counter  flow.  The  average  AT  is  basically  the  same  for  the 
two  flow  types,  but  r\  increases  more  than  linearly  with  AT 
(Fig.  6)  so  we  would  expect  the  parallel  flow  to  offer  more  electric 
power.  We  predict  power  generation  of  33.5  and  32.2  W  for  paral¬ 
lel  and  counter  flow,  respectively,  i.e.  the  flow  type  does  not  have  a 
significant  effect  because  AT  only  varies  in  a  small  range  (where  r\ 
increases  linearly  with  AT).  It  becomes  prominent  in  cases  where 
the  channel  length  is  optimized,  and  has  a  larger  AT.  If  a  channel 
of  optimized  length  can  be  reached,  as  well  as  a  greater  cooling 
along  it,  the  effect  of  flow  types  becomes  more  prominent. 
Fig.  7(b)  captures  AT  between  the  two  fluids  at  the  channel  ends 
as  a  function  of  Lch  which  increases  from  1  to  5.5  m  for  the  parallel 
flow.  Lch  =  5.5  m  is  found  to  be  optimum,  since  AT  approaches  to 
0  K  at  the  reference  conditions  (Tin _coW  =  293  K,  Tin  hot  =  373  K,  and 
AP  =  2  bar)  and  there  is  no  more  available  thermal  energy  to  be 
converted.  Variation  of  power  generation  with  Lch  can  be  seen  in 
Fig.  7(c).  As  expected,  as  LCh  increases,  the  parallel  flow  outper¬ 
forms  the  counter  flow. 

Another  important  factor  is  the  pressure  drop,  AP,  which  influ¬ 
ences  both  the  power  output  and  the  net  power.  The  net  power 


supplied  by  the  TEG  will  be  the  total  power  output  minus  the 
power  consumed  by  the  pump  driving  the  fluid  through  it,  i.e. 

Wnet  =  W'-Wpp  (19a) 

where 

WPP  =  V  AP  (19b) 

is  the  pumping  power,  AP  is  the  pressure  drop,  and  V  is  the  volu¬ 
metric  flow  rate.  It  is  known  that  the  convective  heat  transfer  coef¬ 
ficient  is  changed  by  varying  the  fluid  flow  rate,  which  is  increased 
by  increasing  AP  through  the  channels  of  the  heat  exchanger.  This 
leads  to  a  variation  in  the  thermal  resistances  of  the  heat  exchanger, 
resulting  in  increases  of  AT  and  the  power  output. 

Fig.  8  displays  the  total  power  output  and  the  net  electric  power 
as  a  function  of  AP  with  everything  else  held  constant.  As  AP 
increases,  the  total  power  output  increases  too,  first  relatively 
quickly  and  then  leveling  off.  The  net  electric  power,  on  the  other 
hand,  drops  increasingly.  This  would  seem  to  favor  low  values  of 
AP,  but  not  so  low  that  the  flow  changes  from  turbulent  to  laminar 
(see  below).  Therefore,  AP  along  the  channel  needs  to  be  carefully 
considered  for  the  TEGs  design,  and  is  kept  constant  at  2  bar  in  this 
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Pressure  Drop  (bar) 

Fig.  8.  Variation  of  the  power  generation  and  net  power  with  pressure  drop  at  the 
reference  conditions  ( Tin_Coid  =  293  K,  Tin_hot  =  373  K,  and  parallel  flow). 

study,  which  also  eliminates  issues  with  phase  change  during  the 
heat  transfer  process. 

4.3.  Effect  of  cross-sectional  area  (Ac) 

Next,  variations  in  the  design  of  liquid  channels  and  thermo¬ 
electric  modules  are  assessed  for  fixed,  outside  heat  exchanger 
dimensions  (H  14  mm,  W  12  mm,  L  30  mm).  For  this  purpose,  the 
effects  of  the  cross-sectional  area  of  a  single  channel  Ac  are  evalu¬ 
ated,  assuming  the  same  width  W,  length  L,  and  fixed  operation 
conditions.  Distinction  of  different  TEGs  studied  is  based  solely 
on  the  channel  (thermoelectric  module)  height,  Hch  (HTE),  in  the  z 
direction  (Fig.  3).  In  essence,  the  baseline  TEG  in  Fig.  3  consists  of 
thermoelectric  modules  HTE  =  2  mm  laminated  between  a  hot 
(top)  and  a  cold  (bottom)  channel,  each  of  height  Hch  =  6  mm,  oper¬ 
ating  at  a  specific  total  pressure  drop  AP  =  2  bar.  As  seen  in  Fig.  9, 
the  next  TEG  (Case  2)  consists  of  four  hot  channels  and  four  cold 
channels  assigned  alternately  ( Hch  =  1.5  mm),  with  seven  layers  of 
thermoelectric  modules  ( Hch  =  0.29  mm)  sandwiched  in  between, 
meeting  the  same  pressure  drop  constraint  as  the  baseline  TEG. 
This  iteration  is  repeated  for  sequentially  smaller  distinct  channels, 
as  well  as  thermoelectric  modules,  as  illustrated  in  Table  2.  Eight 
levels  of  mesh  resolution  are  tested  to  check  the  dependence  of 
the  solutions  on  grid  design  for  each  case  in  this  study.  Considering 


Table  2 

Different  TEGs  studied  for  Ac  effect. 


Case 

#  of  channels 

Hch  (mm) 

Hte  (mm) 

1 

2 

6 

2 

2 

8 

1.5 

0.29 

3 

16 

0.75 

0.13 

4 

32 

0.35 

0.06 

5 

64 

0.1875 

0.03 

6 

128 

0.09375 

0.016 

Fig.  10.  AT  across  the  thermoelectric  modules  along  the  flow  direction  are  plotted 
for  the  TEGs  in  the  six  cases. 


the  balance  between  computational  efficiency  and  the  accuracy  of 
results,  suitable  grid  levels  of  resolution  are  chosen. 

Fig.  10  shows  the  variation  of  AT  across  the  thermoelectric  lay¬ 
ers  with  respect  to  the  flow  distance  for  the  different  cases,  indicat¬ 
ing  the  systematic  effects  of  size  as  Ac  moves  from  the  macro-scale 
to  the  micro-scale.  For  Cases  1-3,  AT  remains  significant  until  the 
end  of  the  channel  (30  mm),  but  less  heat  passes  through  without 
contributing  to  power  generation  as  Ac  decreases.  The  flow  dis¬ 
tance  after  which  AT  is  essentially  small,  then  increases  again  from 
an  estimated  1.25  mm  for  Case  4,  through  3.5  mm  for  Case  5  to 
7.5  mm  for  Case  6. 

As  expected,  the  increase  in  effective  heat  transfer  area  at  the 
interfaces  between  channels  and  thermoelectric  layers,  and  the 


Fig.  9.  Isometric  view  (a)  and  inlet  (-x)  view  (b)  of  Case  2  is  shown  in  this  figure  for  Ac  effect. 
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Fig.  11.  uavg  at  outlet  port  as  a  function  of  Ac;  As  Ac  decreases,  the  trendline  with 
three  different  slopes  indicates  three  flow  regions,  i.e.  turbulent,  transition  and 
laminar. 


Case 

Fig.  12.  Total  power  generated  by  the  TEGs  and  Re  with  respect  to  Ac  for  the  TEGs  in 
the  six  cases  at  Tincoid  =  293  K,  TinJlot  =  373  K,  and  AP  =  2  bar. 


reduction  in  individual  layer  thicknesses  as  the  number  of  chan¬ 
nels  and  thermoelectric  layers  increase  tend  to  accelerate  the  heat 
transfer  across  the  TEGs.  However,  the  average  velocity  of  water 
flow  ( Uavg )  in  the  channels  is  lowered,  since  it  is  directly  related 
to  Nu.  This  reduces  the  equilibration  of  temperature  across  the 


Table  3 

Dimension  specifications  for  all  the  five  lengths. 


I  (mm) 

#  of  parallel  unit 

f-chl 

30 

1 

tch2 

7.5 

4 

tch3 

5 

6 

Lch4 

2.5 

12 

tch5 

Case  1 

2.1 

14 

Case  2 

1.6 

19 

Case  3 

1.07 

28 

Case  4 

0.8 

38 

Case  5 

0.5 

57 

Case  6 

0.5 

57 

height  of  a  channel  and  thus  has  a  negative  effect  on  the  heat 
transfer.  Indeed,  as  the  channels  are  getting  sequentially  smaller, 
AT  across  the  individual  thermoelectric  layer  decreases  more 
quickly  along  the  channel;  this  is  not  the  case  in  Fig.  10.  In  fact, 
moving  from  Cases  4  to  6,  uavg  drops  precipitously  as  the  reduction 
in  channel  height  drives  the  water  from  turbulent  flow  to  laminar 
flow  (Fig.  11),  and  AT  remains  high  for  longer  in  spite  of  the  fact 
that  the  effective  heat  transfer  area  keeps  increasing. 

Fig.  12  shows  the  total  power  generation  for  the  six  cases.  The 
corresponding  flow  pattern  is  represented  by  the  Reynolds  number 
(Re)  shown  on  the  secondary  axis  (right).  The  conversion  capability 
is  shown  to  increase  from  Cases  1  to  4  because  of  the  increased 
effective  area  for  heat  transfer  together  with  a  uniform  tempera¬ 
ture  distribution  across  the  height  of  the  channel  caused  by  the 
turbulent  flow.  This  means  that  less  heat  is  lost  by  passing  through 
to  the  end  of  the  channel,  and  a  larger  fraction  of  the  heat  is  trans¬ 
ferred  to  the  thermoelectric  layer  while  AT,  and  thus  r\  is  higher. 
Moving  from  Cases  4  to  6,  however,  the  hydraulic  diameter 
decreases  enough  for  the  flow  to  become  laminar  (reflected  in 
Re)  because  of  the  increased  flow  resistance.  The  conversion  capa¬ 
bility  is  thus  lowered  for  Cases  5  and  6,  even  though  the  effective 
area  keeps  increasing.  We  conclude  that  the  optimal  channel 
height  is  located  at  the  onset  of  the  transition  from  turbulent  to 
laminar  flow,  as  approximated  in  Case  4,  and  determined  by  the 
competition  between  this  (uavg)  and  the  cross  sectional  area  per 
channel,  Ac. 


4.4.  Effect  of  channel  length  (Lch)  &  number 

The  different  effective  distances  along  which  AT  across  the 
thermoelectric  layers  is  still  relatively  large  in  our  calculations 
above  show  that  none  of  the  TEG  designs  considered  are  truly  opti¬ 
mal.  This  can  be  clearly  seen  in  Fig.  10,  particularly  for  Cases  3-6 
where  AT  across  the  TEGs  are  relatively  low  after  x  =  5  mm.  The 


(a) 


(b) 


Fig.  13.  Geometry  specifications  of  Case  1  in  Lch2  at  isometric  view  (left)  and  inlet  (-x)  view  (right).  Hot  and  cold  channels  are  isolated  by  thermoelectric  modules. 
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Fig.  15.  Overall  TEG  efficiency  vs.  Ac  for  LCh\-lCh5  at  rin  co/d  =  293  K,  Tin_hot  =  373  K 
and  AP  =  2  bar. 


six  TEG  designs  are  therefore  considered  again,  this  time  shorten¬ 
ing  Lch  but  keeping  the  total  system  volume  the  same. 

Fig.  13  shows  the  design  from  Fig.  3  where  the  TEG  has  been  cut 
into  Lcji2  =  7.5  mm  long  sections  along  the  length.  The  four  sections 
are  then  placed  next  to  each  other,  providing  for  a  total  width  of 
48  mm.  The  other  five  designs  shortened  (and  widened)  in  the 
same  fashion.  As  indicated  in  Table  3,  the  same  designs  are  also 
cut  into  six  sections  of  Lch 3  =  5  mm  and  12  sections  of  Lch4  =  2.5  mm 
and  combined  into  correspondingly  wider  units  with  more 
channels. 

Fig.  14(a)  shows  the  power  generation  for  each  case  with  the 
different  lengths  and  corresponding  numbers  of  parallel  channels. 
Shorter  lengths  and  more  parallel  channels  provide  for  a  reduced 
flow  resistance  and  a  greater  effective  area  for  heat  transfer  near 
the  inlet  where  the  AT,  and  thus  the  thermoelectric  conversion 
efficiency,  is  largest.  For  a  given  length,  the  total  power  goes 
through  a  maximum  value  determined  by  the  competition 
between  heat  transfer  area  and  flow  nature  as  before.  The  reduced 
flow  resistance  is,  however,  also  accompanied  by  an  increase  in  the 
Reynolds  number,  i.e.  there  should  be  a  tendency  for  the  maximum 
to  shift  towards  smaller  channel  heights.  One  more  set  of  calcula¬ 
tions  is  therefore  performed  in  which  the  optimum  channel  length 


is  identified  for  each  of  the  six  cases  (Table  3).  The  results,  included 
in  Fig.  14(a),  show  a  shift  of  the  optimum  to  near  Case  5  with  an 
optimized  Lch  as  the  corresponding  Re  increases  from  422  to 
1191  to  3533  (Fig.  14(b)),  and  the  flow  thus  changes  from  laminar 
to  turbulent,  for  this  channel  height  (0.1875  mm).  This  should  be 
near  the  true  optimum  design  for  the  present  operating  conditions. 

4.5.  Overall  TEG  efficiency  (rj) 

Fig.  15  shows  the  overall  TEG  efficiency  ri  (the  portion  of  energy 
that  can  be  extracted  from  the  liquid)  with  respect  to  Ac  for  all  the 
lengths  studied.  Compared  to  Fig.  14(a),  the  more  power  genera¬ 
tion  associated  with  optimizing  the  TEG  configuration,  the  higher 
r\  can  be  achieved.  Case  5  in  LChs  shows  the  highest  r\  of  up  to 
3.5%.  This  should,  however,  only  be  considered  as  an  ideal  value, 
since  several  additional  considerations  need  to  be  incorporated 
into  the  proposed  thermoelectric  model  for  further  studies. 

5.  Conclusion 

Two  thermoelectric  models  are  applied  to  assess  the  device  per¬ 
formance  of  the  TEGs,  and  the  results  show  that  the  Coupled-field 
model  performs  better  compared  to  the  Simple  model  when  vali¬ 
dated  using  prior  experiments  [13].  The  Coupled-field  model  is 
capable  of  representing  the  energy  balance  at  the  interface  (TE  ef¬ 
fects)  and  the  predicted  power  generation  is  within  8%  error  com¬ 
pared  to  the  experiment  validation.  This  is  a  much  more  improved 
result  compared  to  most  thermoelectric  simulations  available  in 
the  literature. 

The  simulations  based  on  the  Coupled-field  model  offers  signif¬ 
icant  insight  into  the  systematic  effects  of  fluid  heat  exchanger  de¬ 
sign  parameters  on  the  efficiency  of  a  generic  thermoelectric 
generator  (TEG).  Not  surprisingly,  increasing  the  inlet  temperature 
(Tin_ hot)  of  the  hot  liquid  leads  to  increases  not  only  in  the  energy 
input,  but  also  in  the  conversion  efficiency  with  which  this  can 
be  converted  to  electricity.  Parallel  flows  of  the  hot  and  cold  liquids 
lead  to  greater  overall  TEG  efficiency  as  compared  to  counter  flows. 
The  pumping  power  required  to  maintain  a  pressure  drop  (AP) 
along  the  channels  is  significant,  more  than  increasing  the  conver¬ 
sion  capability,  but  can  counteract  the  effect  of  a  higher  pressure 
drop  on  the  net  power. 

In  general,  it  is  critical  to  maintain  turbulent  flow  in  channels.  A 
larger  number  of  channels  at  the  micro  scale  provide  for  more 
interfaces  with  thermoelectric  layers  between  them  and  smaller 
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temperature  gradients  across  the  channels,  but  eventually  the  in¬ 
creased  flow  resistance  leads  to  a  transition  to  laminar  flow  and 
lowers  the  power  output,  unless  the  pressure  drop  is  raised.  The 
power  generation  is  therefore  optimized  by  simultaneously  mini¬ 
mizing  channel  dimensions  (height  and  length)  as  much  as  possi¬ 
ble  without  allowing  (i)  the  flow  to  become  laminar  or  (ii) 
significant  thermal  energy  to  exit  out  the  other  end  of  the 
channels. 

It  follows  that  further  optimization  may  be  achieved  by  intro¬ 
ducing  a  stirred  flow  to  mix  the  stream  at  low  Re.  Future  work  will 
also  incorporate  enhanced  mixing  using  passive  microstructures 
and  high  thermal  conductivity  nanofluids  into  the  numerical 
model. 
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